Statistics of resonance states in open chaotic systems: A perturbative approach 
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We investigate the statistical properties of the complexness parameter which characterizes uniquely complex- 
ness (nonorthogonality) of resonance eigenstates of open chaotic systems. Specifying to the regime of weakly 
overlapping resonances, we apply the random matrix theory to the effective Hamiltonian formalism and derive 
analytically the probability distribution of the complexness parameter for two statistical ensembles describing 
the systems invariant under time reversal. For those with rigid spectra, we consider a Hamiltonian character- 
ized by a picket-fence spectrum without spectral fluctuations. Then, in the more realistic case of a Hamiltonian 
described by the Gaussian orthogonal ensemble, we reveal and discuss the role of spectral fluctuations. 



PACS numbers: 05.45.Mt, 03.65.Nk, 05.60.Gg 
I. INTRODUCTION 

In the domain of wave or quantum chaos IH, open systems 
are currently actively investigated both from experimental and 
theoretical points of view (see Refs. i, H for recent re- 
views). Openness may be due to various physical mechanisms 
such as bulk absorption, coupling to the environment through 
physical channels as well as dissipative or radiative bound- 
ary conditions. Whatever the mechanism, openness results in 
spectral broadening ranging from the perturbative regime of 
non-overlapping (isolated) resonances to the so-called Erics- 
son regime of strong overlap. These mechanisms and their 
related spectral effects have been experimentally studied in 
various context: in microwave cavities flsSSllllijn opti- 
cal microcavities 11 [Bin, and in elastodynamics 11121 11311 . 

The most salient feature of open systems is the set of res- 
onances which are quasibound states embedded in the con- 
tinuum. A natural way to address them analytically is via 
the energy-dependent scattering matrix, S{E). Following the 
Heidelberg approach jT4|, the poles (i.e., resonances) of the S 
matrix turn out to be the complex eigenvalues of an effective 
non-Hermitian Hamiltonian Jics, whereas the bi-orthogonal 
eigenvectors of the latter determine the corresponding reso- 
nance states (quasimodes). Universal properties of resonance 
scattering in the chaotic regime can then be analyzed by ap- 
plying random matrix theory (RMT) that amounts to replacing 
the actual non-Hermitian Hamiltonian with an RMT ensemble 
of the appropriate symmetry class f\3\. The main advantage 
of such an approach is that it treats on equal footing both the 
spectral and scattering characteristics of open chaotic systems 
as well as that it is flexible enough to incorporate other imper- 
fections of the system, e.g., disorder and losses [jst]. 

By now, complex eigenvalues of such non-Hermitian ran- 
dom matrices have been studied quite systematically lfl6l [itI 
[isll . However, the statistical properties of the correspond- 
ing (left and right) eigenvectors are less understood. Quite 
a substantial progress in this direction has been achieved by 
Schomerus et al. lfl9ll . who studied mainly the systems with 
broken time-reversal symmetry. Other analytical results for a 
few physically interesting particular cases have also been re- 



ported in the literature recently ll20l l2ll l22l l23ll . Components 
of eigenvectors appear as residues of the S matrix at resonance 
positions and the understanding of their properties is thus im- 
portant for many applications. For example, nonorthogonal- 
ity of resonance eigenstates yields the enhancement (the so- 
called Petermann factor) of the line width of a lasing mode 
in open resonators |[l9[| and influences branching ratios of nu- 
clear cross-sections 11241 l25ll . It features also in the particle 
escape from the scattering region ll26ll as well as in dissipative 
quantum chaotic maps luTIl . 

This paper focuses on spectral and eigenvector statistics 
of such non-Hermitian random matrices describing effective 
Hamiltonians of open chaotic wave systems whose closed 
limit displays time-reversal symmetry (TRS). In this case, the 
quasimodes correspond to the complex-valued eigenvectors 
of Tieff- To characterize this complexness, it is convenient 
to introduce lfl2i l28ll the ratio of the variances of the imag- 
inary and real parts of the eigenvector as a single statistical 
parameter, hereafter called the complexness parameter ll23ll . 
One should note that this parameter is characteristic of the 
degree of non-orthogonality of the complex modes and, there- 
fore, is closely connected to the Petermann factor mentioned 
above lfl9ll . Other studies have considered the phase rigidity, 
another related parameter, introduced to characterize the de- 
gree to which a general scattering wave function is complex 
||22[ I29I1 . Both parameters are straightforwardly deduced from 
one another when the phase rigidity is calculated for a sin- 
gle eigenvector. The main advantage of considering the com- 
plexness parameter is to reveal a physical connection between 
spatial and spectral statistics ll5i l23ll . 

In what follows, we study the probability distribution of 
the complexness parameter for a generic weakly open chaotic 
system and its connection with the distribution of resonance 
widths within the RMT approach. At the first stage, we derive 
an expression for the complexness parameter in the weak cou- 
pling regime and establish a general relation between its aver- 
age and width fluctuations. Then accounting for the essential 
statistical feature of spectra in chaotic systems, namely, spec- 
tral rigidity, we investigate the case of a system whose closed 
limit is described by a pure picket-fence spectrum. An ex- 
act analytic prediction for the probability distribution of the 
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complexness parameter is derived, depending on only two 
parameters; the number of open scattering channels and the 
mean resonance width. Finally, we consider the more realistic 
case of systems modeled by the Gaussian orthogonal ensem- 
ble (GOE). We derive an analytic expression for the probabil- 
ity distribution of the complexness parameter in this case and 
discuss the effect of spectral fluctuations. 



n. EFFECTIVE HAMILTONIAN FORMALISM 

A. Scattering approach 

Open wave systems are commonly described using the so- 
called projection formalism Izsllsoll . The exterior coupling is 
modeled by M scattering channels connected to N levels of a 
closed system. The coupling to the environment turns modes, 
with a infinite life time, into resonances, with a finite life time. 
Being initially introduced in nuclear physics, this formalism 
has been later applied successfully to wave billiards jsill for 
which antennas and absorption can be described by scattering 
channels ll32ll . In this approach, the resonance part of the 5- 
matrix is given by: 



S{E) 



1 



E — Hcs 



-V, 



(1) 



where V is the coupling matrix of size N x M, the elements 
of this matrix couple the n-th level to the c-th scattering 
channel. The poles of S are given by the eigenvalues of HcS- 
Assuming an independence of the coupUng elements from 
the energy and neglecting direct processes IIMl . the effective 
Hamiltonian of the open systems is represented as follows: 



(2) 



where H is the Hamiltonian of the closed system and the anti- 
Hermitian part ^VV"^ describes coupling to the channels. In 
the case of the systems with preserved TRS considered below, 

is a real symmetric matrix of size N x N and V is also real. 
As usual, the limit ^ oo is to be finally taken. 

Since Hcff is a non-Hermitian operator, the eigenvalue 
problems Hcffli^n) = ■S'rjV'n) and (^AnlT^cff = Sn{ijJn\ de- 
fine two sets of a priori independent eigenvectors, called right 
{ I f/'ri ) } and left { ( V'n | } eigenvectors associated to the same set 
of eigenvalues {£„}. These eigenvectors form a bi-orthogonal 
set which satisfies conditions of orthogonality, (f/'nl^m) = 
Snm, and completeness, lV'ri)(V'Ti| = 1- Making use of 
the right eigenvectors, the diagonalization of HcS then reads: 



HofflV'n) = (En - -r„)|V'„) 



(3) 



where En and r„ are, respectively, the energy and the width 
of the n-th resonance. Due to TRS present, HeS is a complex 
symmetric matrix; hence, the left and right eigenvectors are 
related by the transpose, = (IV'n))^ Hit]. 

The coupling to continuum, as described by the imaginary 
part of Hcff, turns real eigenfunctions of the closed system 



into complex quasimodes of its open counterpart. In order 
to measure their complexness, we define the complexness pa- 
rameter Qn as follows: 



(4) 



where i/j^ is the i-th component of the eigenvector (we note 
that the complexness parameter can be equivalently defined by 
means of the left eigenvectors). It is worth noting here that in 
contrast to the related Petermann factor lfl9ll . which is defined 
for a fixed value of the given resonance width, no additional 
constraints are imposed on (|4]i. In chaotic systems, reveals 
strong mode-to-mode fluctuations, which we describe through 
its probability distribution function to be derived below. 



B. Statistical assumptions 

Within the RMT approach, the universal statistical prop- 
erties of closed chaotic systems with preserved TRS are de- 
scribed by GOE ill]. In this ensemble the joint probability 
distribution, P{{Ei}), of the levels (the eigenvalues of H) is 
induced by a Gaussian distribution of the random real sym- 
metric H with zero mean. The exact expression for P{{Ei}) 
is well known to have the following form: 



Pi{E^})<X n \En-E,n\cxp 
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(5) 



Here, we have chosen the variance of H such that it yields the 
mean level spacing A = 1/iV at the spectrum center, E = {). 

The energy levels, as defined by Eq. (|5]l, exhibit a lin- 
ear level repulsion. As a result, the energy spectrum dis- 
plays spectral rigidity which restrains the spectral fluctuations 
around the mean. This important feature can approximately be 
taken into account within the so-called picket-fence model of 
equidistantly spaced levels |34j]. The usefulness of this model 
is in its simplicity that allows one to treat various resonance 
phenomena analytically, see, e.g., Refs. Issl l36l IjtIi . Here, 
we employ this model to single out a contribution to q^-^ due to 
fluctuations of the resonance widths. 

As concerns the coupling amplitudes, the results are known 
to be model independent on statistical assumptions on as 
long as the number o f op en channels is small compared to 
that of the levels Issl l39ll . The coupling amplitudes may be 
equivalently chosen as fixed ifTill or random Issll . In order 
to preserve orthogonal invariance of TicS under (complex) or- 
thogonal transformations lO , we consider the l^'s as real 
Gaussian random variables with zero mean and 



i2KA/TT)Sn 



2 



(6) 



Henceforth, (• • •) stands for the statistical averaging over the 
ensembles. The coupling constant k determines a transmis- 
sion coefficient T = 1-| (S*) p = 4k/(1+k)2 of the channels 
(assumed to be statistically equivalent). The cases of T <C 1 
or T = 1 correspond, respectively, to weak or perfect cou- 
pling. In the weak coupling regime considered below, k ^ 1, 
all the resonances are almost isolated and (F) <C A. 
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III. PERTURBATIVE APPROACH 
A. Complexness parameter in the weak coupling regime 

We now derive an expression for the complexness param- 
eter of the eigenvectors for weakly overlapping resonances. 
The matrix representation of Hcs in an arbitrary basis {\n)} 
of the Hilbert space spanned by eigenvectors of H reads: 



Ti-cs — 



N 

E 

n,p— 1 



TV M 



n,p—l c—1 



As we focus on the weak coupling regime, the imaginary part 
may be viewed as a perturbation of the Hamiltonian of the 
closed system. The repulsion of the energy levels exhibited 
by the systems under investigation allows us to consider the 
eigenenergies of H as nondegenerate. One can therefore ap- 
ply first-order perturbation theory to obtain from (|3]l the eigen- 
values and the eigenvectors of Hcs straightforwardly. The 
eigenvalues read En — f r„, where the E^s are the eigenval- 
ues of H and the widths r„ are given by: 



M 

E 



The perturbed eigenvectors of TCcS written in the eigenbasis 
{|0„)} of H are easily found as follows: 



\lpn) = -«E 



2{En — Ep) 



(9) 



Splitting then the real and imaginary parts of \tpn), the com- 
plexness parameter (|4li of a given eigenvector reads 



E 



r2 

np 



4(S„ - Ep 



where we have introduced r„p = X^c^i ^n^p- These quan- 
tities are responsible for the coupling and interference of the 
resonance states due to the common decay channels Issll . 

In what follows, we study the statistical properties of the 
complexness parameter ( fTOb for H being described by a 
picket-fence or belonging to GOE. It is worth noting here 
that expression (fTOl i is a sum of correlated random variables 
which, therefore, does not obey the standard central limit the- 
orem. Statistics of a similar kind of objects appears, e.g., in 
the study of the parametric level dynamics ("curvature") ll40ll 
and in the context of interference effects in neutron scattering 
from compound nucleus Ii4li1 . 



(10) 



B. Rescaled parameters and their statistics 

The complexness factor ( fTOb contains two contributions of 
distinct types, one is due to the internal levels and the other 
is due to the coupling matrix elements T^p- From a statis- 



one another We note, however, that the levels En are mutu- 
ally correlated. The quantities r„p's, unlike the original am- 
plitudes Vn, are also not statistically independent. Although 
their joint distribution can be found from (|6]l, the resulting ex- 
pression is quite complicated ll42ll . being of little practical use 
for actual calculations in the present context. 

To overcome the difficulty of averaging over the coupling 
amplitudes, we follow Sokolov and Zelevinsky jssll and treat 
an arbitrary matrix element r„p as a scalar products between 
A/-dimensional vectors V„ and Vp of the coupling ampli- 
tudes {Vf } associated with the levels I = n and / = p. This 
suggests a natural parametrization for r„p in terms of the an- 
gles 9np between the pairs of these N vectors. 



Tnp = (V„ • Vp) = v/raVcos0„p. (11) 

The main advantage of this representation is that the angles 
6np are mutually independent and also independent of r„. 
The probability distribution of any angle (for M > 2) can 
be easily found to be given by the expression for a solid angle 
in an M-dimensional space Il33ll : 

r(7\//2) 



PMiO) 



V^r((M-i)/2) 



(12) 



(8) Note that r„p = yT\JT^ at AI = 1. As concerns the lengths 
of these vectors, i.e., the widths dH), these are well-known to 
be independent and distributed according to 

1 



PAlil) 



2M/2p(^^/2) 



^A//2-lg-7/2 



(13) 



Henceforth 7„ = r„/(T^ stands for the dimensionless widths. 
This distribution function has the mean value (7) — M and 
the variance 

var(7) = 2M = A . (14) 

Thus the widths cease to fluctuate as the number of open chan- 
nels grows, with the average width being kept fixed. 

It is now convenient to express all the quantities in their 
natural units and to consider a rescaled complexness parame- 
ter Xn defined as follows: 



_ A2 2 

Xn = —Tin 



In 



pi^n 



Epf 



(15) 



(16) 



where we have introduced the following quantities 

Zp = 7p COS^ Onp ■ 

Zp may be given a geometrical interpretation as (a square of) 
the projection of the vector cr~^Vp along the direction given 
by the vector V„. These projections are statistically indepen- 
dent, as is obvious from the above discussion. The probability 
distribution of any projection follows readily from Eqs. ( fT2b 
and ( HjI ). Performing an integration first over 7 and then over 
9 in the definition P{Z) = {5{Z - 7003^ 9)), one finds 



P{Z)^ 



1 



-Z/2 



(17) 



Thus, surprisingly, the distribution of Zp is independent of M, 



tical point of view, these two are statistically independent of being given by the Porter- Thomas law at any M > 1. 
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C. Average of X and width fluctuations 

A general expression of the average value of the complex- 
ness parameter X can be readily found from Eq. ( fTSl l by mak- 
ing use of the mutual statistical independence between the 
widths {7„}, the projections {Zn} and the levels {En}. Not- 
ing that (7) = M and {Z) = 1 = (7) ^cos^ 9^, one obtains 



{X) = Alf , 



(18) 



where the factor / depends on the statistical properties of the 
energies of the closed system only, 



A2 



A{En - Ep) 



(19) 



It is important to note that, generally, the nonzero values 
of the complexness parameter are solely due to fluctuations 
of the resonance widths. Indeed, in the extreme case of all 
the widths being equal, the anti-Hermitian part of Hcs gets 
proportional to the unit matrix and, as a result, the complex 
(biorthogonal) eigenvectors become essentially real 1I23I1 . It is, 
therefore, instructive to take this explicitly into account and, 
in view of relation (fT4l i. bring Eq. (fTsT l to the form: 



{X) = (var(7) . 



(20) 



This expression relates the average complexness parameter to 
the natural measure of the width fluctuations, its variance. 

Strong correlations between the complexness parameter 
and the spectral widths are already known. The proportion- 
ality between ^ {X) and the average value of the fluctuat- 
ing part of damping was recently found experimentally in a 
chaotic microwave billiard at room temperature, where this 
was also explained heuristically using a ray picture based on 
the ergodic character of the wave system 15|l. Then this pro- 
portionality was established in ll23ll using a two-level RMT 
model and considering M 3> 1 that was relevant for this ex- 
periment. Expression (|20] | readily provides this feature, in 
view of 



, {X) = (7) yfJjM, at any N and M. On the other 
side, it captures fluctuation properties of the widths properly, 
e.g., yielding the vanishing (X) in the absorptive limit of 
many weakly coupled channels with the average total width 
kept fixed, due to the vanishing variance (fT4l i. Therefore, 
we believe that relation (|20| | is a general feature of weakly 
open chaotic systems with non-degenerate spectrum in the 
perturbative regime. Figure 1 supports this suggestion through 
numerical simulations of the picket-fence and GOE models 
(with the details being given later in the next section). 

A remark on the proportionality factor / is appropriate here. 
In the RMT limit N 00, this factor may be represented 
as follows / = ^ J^ds s^^R2{s), where i?2(s) is the two- 
point correlation function of the RMT. The main problem of 
the GOE case, already mentioned in lfl9ll23ll . is an 'infrared' 
logarithmical divergency of / due to i?2(s) ^ s at s — > 0. 
Practically, this divergence can be regularized by introducing 
a cut-off at small s, s > e, see Fig. 1. Without this cut-off the 
expression of the complexness parameter obtained using first 
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FIG. 1: The average rescaled complexness parameter versus the 
width variance for the GOE and picket-fence models. The sym- 
bols correspond to the results of numerical simulations performed 
at Af = 1, 2, 3, 5, 10 (see the text for details). The linear depen- 
dence predicted by Eq. i20i is represented by the solid line. In the 
GOE case (a), the proportionality factor is given by the regularized 
expression /e = ^ ST^ ds s^^ R2{s). The results obtained with three 
different values of the cut-off e are shown. In the picket-fence case 
(b), / = 7r^/12 as exactly given by Eq. l |23t . 



order perturbation theory ( fTSl l does not yield finite moments, 
thus demanding for the characterization of fluctuations of X 
by means of its probability distribution. 



IV. DISTRIBUTION FUNCTION 

The probability distribution function of the rescaled com- 
plexness parameter Xn is defined as follows 



VAiiX) = {S{X-Xn)) 



(21) 



where the statistical averaging over the levels, the widths and 
the projections is performed with the help of Eqs. (|5]), (O, 
and ([TtT i. respectively. In the weak coupling regime, function 
( l2n i depends only on the number AI of open channels. 

It is instructive first to consider the case of the completely 
rigid spectrum, which may be viewed as an approximation of 
the GOE spectrum where the fluctuations are neglected. 
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A. The picket-fence model 

In this model the eigenenergies of the closed system are 
equally spaced, i.e. En — En±k = i/cA, and the eigenvector 
components are random Gaussian variables. The complexness 
parameter is then given by 



fc#0 



(22) 



This expression does not have any divergence problems of the 
GOE case, thus statistics of (l22l i can be also characterized by 
its moments. In particular, the average value is easily found to 
be exactly given by Eq. ( |20] |. with the factor / being 



f=y- 

J Ah 



TT 

12 



(23) 



Figure 1(b) illustrates the dependence (X) ~ ^var(7). 

We now derive an exact expression for the probability dis- 
tribution V^iX) in the picket-fence case. First we substitute 
in the definition ( ISTT i the Fourier representation of the delta 
function, S{X - X„) = / l^e*'^^'^-^"), where X„ is given 
by Eq. (l22T i. Then the integration over the projections with 
the help of Eq. (T7\ becomes trivial, yielding 



dto 
2^ 



M 



(7)nr 



(24) 



fc=i 



The infinite product here can be evaluated explicitly II43I1 . 
Making use of the explicit expression ( fT3] l for Pa/ (7) and ap- 
plying the change of variables 7 = 2|2p, Eq. ( |24| ) can then be 
cast in the following form; 



1 



dhJ 



V{M/2) J 27r 



— 1 1 1 1 1 1 r- 

10-1 



10- 



2 3 4 5 6 7 



lllTTTTTrrr j ^ 




M-l AujX- 



,2 yJlOJZ-K 



sinh(vzwz7r) 



(25) 



As one can easily check, this expression is properly nor- 
malized to unity. It is also worth noting that the integrand of 
Eq. ( I25] ) is an analytic function in lo except for the poles lo- 
cated on the upper part of the imaginary axis at Uk ~ i{k/ zY, 
k = 1, 2, . . . , 00. This readily implies that V\l{X) = at 
X < identically. 

The details of the subsequent calculations of V^jiX) are 
given in AppendixlAl The final expression reads: 



Vt,{X) = 



2^{VxY'^/ 
r(M/2) 



2-1 



dz . . Jm/2-i(2VXz), (26) 
smn(z7rj 



FIG. 2: (Color online) The distribution of the rescaled complexness 
parameter for the picket-fence model at M = 1, 5 and 10 (top, mid- 
dle and bottom, respectively). The analytical result l i25b is plotted in 
the solid line while the histograms correspond to numerics. Insets 
show the tail of the distribution in a semi-log scale. 



with (x) being the Bessel function of order ly. In the case of 
an odd number of channels, M ^ 2n + 1, n = 0,1, . . . , this 
expression can be integrated further to yield an attractively 
simple formula 



1 



2r(n + l/2) \ dXj cosh2(Vx; 
In particular, the single-channel distribution Vi^X) reads 



(27) 



Vf{X) 



1 1 



2V^cosh^(V^) 



(28) 



6 



It is interesting now to study in details the case of the large 
number of weakly open channels M ^ 1. In view of the 
scaling ( fTsl l, we consider the limiting probability distribution 
of .T = X/M defined as 



p{x) = lim MVm{Mx) . 



M 



(29) 



Expression (126b is actually not very convenient for evaluating 
this function. However, one can note that in the limit con- 
sidered, the distribution Pj^jIj), Eq. ( fTlT l, tends to the Dirac 
distribution, 5(7 — M). Then, starting from Eq. ( |24] |. the in- 
tegration over 7 is trivial and the probability distribution of x 
reads: 



1 



(30) 



Using the residue theorem, one readily gets: 



Ppf(x)=4^(-1)^ 
fe=i 



(31) 



and finally 



ax 



(32) 



where -d^ is a Jacobi theta function Il43ll . 

The above analytical predictions concerning the average 
value of the complexness factor and its probability distribu- 
tion have been checked through numerical simulations of ran- 
dom matrices, see Figs. 1 and 2. Numerical simulations are 
based on the diagonalization of the effective Hamiltonian (|2]i 
viewed as a random non-Hermitian matrix. We have consid- 
ered resonances in the bulk only, i.e. resonances with a large 
number of neighbors on the left and on the right of the spec- 
trum. This restriction is introduced to neglect the edge effects 
whose contribution tends to vanish as ^ 00. 

The picket-fence Hamiltonian is built such that the eigenen- 
ergies are equally spaced and the eigenvectors are random 
Gaussian variables. This is readily done by following a proce- 
dure adapted from ll44ll where the authors used it to generate 
the POE ensemble. Thus, in a basis deduced from its eigen- 
basis through an arbitrary orthogonal transformation O with 
random Gaussian variables, the Hamiltonian H is given by: 



H ^Odmg{En}0' 
where En = n/N, such that A = l/N, and 



(33) 



(34) 



Statistics were performed with 100 matrices of size 
lOOOx 1000. In order to make the calculated distributions in- 
sensitive to edge effects, 100 levels at each end of the spec- 
trum were discarded. In all the simulations the mean spectral 
width is kept fixed and equal to (F) / A = 10^^. 
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FIG. 3: (Color online) The distribution of the rescaled complexness 
parameter for the GOE model at A/ — 1,5 and 10 (top, middle and 
bottom, respectively). The analytical result l l35t is shown in the solid 
line and compared to that l |25t of the picket-fence case (dashed line) 
while the histograms correspond to numerical simulations. Insets 
show the tail of the distribution in a log-log scale. 



B. The GOE model 

The probability distribution in the GOE case can be found 
by makinguse of group integration methods and results ob- 
tained in 111911 . Outlining the details of the computation in Ap- 
pendix|B] we state the final result here: 



A/ 



TT^M l + 7r2(3 + M)/(4X) 
24X2 [l + 7r2/(4X)]^V2+2 ■ 



(35) 



To check our findings, the same kind of numerical simula- 
tions as in the picket-fence model have been performed. The 
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closed Hamiltonian H now belongs to GOE, its elements be- 
ing defined by their first two moments; 

where N is the size of the matrix. Like in the picket-fence 
case, the normalization is chosen such that A = Statis- 
tics were obtained with 150 matrices of size 1000x1000. 
Only levels near E =^ for which spacings deviate less than 
5% from A were kept. The agreement between numerical and 
analytical results is flawless, as shown in Fig. 3. 

The comparison between the probability distribution of X 
in the picket-fence model and for GOE illustrates the effects of 
the fluctuations of the spectrum on the complexness parame- 
ter The maximum of both distributions are close to each other 
This is mainly due to the spectrum rigidity in both ensembles. 
But at large X the statistical weight is larger for GOE than 
for the picket-fence model. This difference is introduced by 
the behavior of the levels at small distance: the spacing of 
two eigenenergies can be very small, the corresponding con- 
tribution to the complexness parameter is large, then the tail 
of Vm{X) is larger for GOE than for the picket-fence model. 
This feature is most explicitly seen by comparing the corre- 
sponding limiting distributions at M ^ 1. The distribution 
( |29] l of X is easily obtained from Eq. ( [35] ) and reads 

2 2 / 2 \ 

In contrast to the asymptotic exponential behavior in the 
picket-fence case, Ppf (x) cx e~^^', see Eq. ( [3T1 ), the tail of the 
distribution (|37] | follows a power-law decay: pgoo(a;) oc 

V. CONCLUSION 



may be considered as a sensitive probe of the crossover from 
localized to extended states in open disordered systems |47]. 
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APPENDIX A: DERIVATION OF EQS. <26) AND 0?) 

We first note that the integrand of Eq. dZST l is a symmet- 
ric function in z that allows us to restrict the z-integration to 
the positive axis. Then we deform the contour of integration 
over u) from the real to imaginary axis by putting O — ius. 
Performing after that the scaling transformations of the inte- 
gration variables, first z z/\/f7 and then Q fl/X, and 
interchanging the order of integrations over z and ft, we may 
cast Eq. ( |25] l in the following form 

' r(A//2) J sinh(7rz) 



-\-ioo 

X / ^f7-M/2gO-x.Vn (Al) 
J 2m 

— zoo 



In this paper, we have studied the statistics of complex 
wavefunctions associated to the resonances of weakly opened 
wave chaotic systems with the preserved time-reversal sym- 
metry. More specifically, in the perturbative regime, we have 
considered the case of the completely rigid spectra defined 
through the picket-fence model and that of the GOE display- 
ing spectral fluctuations. One of the key features of this study 
relies on the proportionality between the average of the com- 
plexness parameter and the variance of the resonance widths, 
which we believe is valid for generic nondegenerate spectra. 
We have also derived the exact probability distribution of the 
complexness parameter in these two cases. 

To check the validity of the present results, recent exper- 
iments in elastodynamics are available. In particular, in the 
case of vibrating plates, a complete knowledge of the eigen- 
functions can be obtained through noninvasive measurements 
||4^ even for moderate overlap of resonances. Indeed, the 
understanding of the statistics of eigenfunctions beyond the 
perturbative regime still remains an open problem. (We note 
that some relevant interesting numerical results for microwave 
billiards with large openings were recently reported ll46ll .) Fi- 
nally, one should also note that the complexness parameter 



To calculate here the last integral over 57, we expand e ^^^/^^ 
into a series and evaluate the result termwise 



E 



-\-ioo 

k\ I 2m 



i-Xz^) 



i-WXz) 



\2]k 



(A2) 



where we have used f^fl-^e^^ = l/^M- Making 

now use of the well-known series representation for the Bessel 
function ll43ll . one can immediately recognize the r.h.s. of 
dAa to be equal to Jm/2-i(2\/Xz). Collect- 

ing all the factors together, we finally arrive at Eq. ( |26] l. 

Further progress is possible in the case of odd M . It is 
instructive first to start with the case of M = 1, which turns 
out to play the central role in this calculation. We may use 
the known relation J_i/2(z) = ^J2z / tt cos{z) / z in this case 
m, thus (V^z)i/2j_i/2(2V^z) = ■^cos(2VXz), that 



8 



allows us to perform the integration in Eq. ( lAll ) analytically: 



dz 



z f dQ 



sinh(z7r) J 27ri 



= d 



zcos{2VXz) 1 1 
' sinh(z7r)0F ~ 4V^cosh^(V^) 



(A3) 



Taking now into account the (omitted) factor 2y^7:/X, we ob- 
tain vf{X),Eq. (|28l). 

The general case of odd A/ = 2n + 1 may be reduced 
to that of M = 1 considered above, if one notices that the 
term z^^/fi^^/^e^''^^'/^ in the integrand of Eq. dAll ) can be 
generated by a differentiation with respect to X as follows: 



^-xz^/n _ 



d 



dX / 



(A4) 

Substituting this representation into Eq. (lAlb and changing 
the order of the integrations and differentiation there, we see 
that the resulting integral is already given by Eq. (|A3t that 
readily yields the expression dZTj ) of Sec. llV A1 



and {En} are taken from the GOE. They found the following 
expression for P{A, B): 



We note that the above expression ( |B3t was obtained in 
ifioll for the particular case of one open channel. The key fact 
which allows us to apply this result to our ^/-channel case is 
the representation ( fTSl ) in terms of projections with the distri- 
bution ( [TtI i. The later corresponds to the Gaussian distribution 
(IB 2b with K = 1 and -^oip = Zp, thus giving a connection 

X = ^jB. Correspondingly, the distribution function of X 
in the GOE case can be found from 



rfr{x) = (5ix^^jB) 



(B4) 

by averaging over A, B and 7. Substituting the explicit form 
(|B3b , it is convenient first to integrate out B that yields 



/27r 
"12" 



d7^M(7)(a7)'/'e-'^'^/2 



APPENDIX B: DERIVATION OF EQ. 

We use the recent result by Schomerus et al. ifTgll . who 
calculated the joint probability distribution P{A, B) of 



Ep — En 



p^n 



{Ep~En?' 



(Bl) 



where {ctp} are the statistically independent real Gaussian 
variables distributed according to 



2KAa2 



-7r2a^/(2,cA) 



(B2) 



j dA{l + 7r^A^)e-^-^/^^^"''\ (B5) 



with a = 7r^/4X. The Gaussian integration over A is now 
straightforward and gives 



^M°(^) - 24r(M/2) X2 



where we have substituted expression ( fTST l for Pm{j)- The 
remaining integration yields Eq. ( |35] |. 
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